Microscopic origin of the linear temperature increase of the magnetic susceptibility of 

BaFe 2 As 2 
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Employing a combination of ab initio band structure theory and dynamical mean-field theory we 
explain the experimentally observed linear temperature increase of the magnetic susceptibility of 
the iron pnictide material BaFe2As2. The microscopic origin of this anomalous behaviour is traced 
to a sharp peak in the spectral function located approximately 100 meV below the Fermi level. This 
peak is due to the weak dispersion of two-dimensional bands associated with the layered crystal 
structure of pnictides. 

PACS numbers: 74.70.Xa, 71.27. +a, 71.10.-w 



I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity in the iron pnictides 1 - the unusual electronic prop- 
erties of this class of materials has attracted considerable 
attention^ ^ 11 ' 14 ' 22 ! 23 . The interest was stimulated by 
the fact that the new superconductors share several sim- 
ilarities with the well-studied, but still not well under- 
stood, high-T c cuprates. Firstly, both classes of super- 
conductors crystallize into a layered structure. Secondly, 
in most cases the parent compounds of the pnictides 
are not superconducting, and superconductivity emerges 
only under doping or pressure and is associated with 
the suppression of antiferromagnetic (AFM) order. How- 
ever, in contrast to the cuprates whose parent compounds 
are Mott insulators, the parent compounds of the pnic- 
tides are multiband metals. The nature of the magnetic 
ground state of the parent compounds is also different: 
in the cuprates it corresponds to a Neel-type order of 
a Mott-Hubbard insulator, while in the pnictides mag- 
netism is associated with a nesting-induced spin density 
wave^ (SDW). 

The magnetic properties of pnictide materials show 
anomalous behavior even in the paramagnetic state. An 
unusual linear temperature increase of the uniform mag- 
netic susceptibility was reported in the parent compound 
BaFe2 As^ as well as in stoichiometric and fluorine-doped 
LaFeAsO 6 . It is now well established that the linear in- 
crease of the magnetic susceptibility with temperature 
is a general property of all pnictide superconductors for 
temperatures above the SDW transition. Nevertheless, 
no consensus has been reached so far about the origin 
of this phenomenon. To date several mechanisms were 
proposed to explain the observed T-dependence in the 
pnictides. Wang et ai~ and Zhang et alJ- suggested that 
the linear- T behavior is a consequence of strong antiferro- 
magnetic fluctuations present above the SDW transition 
temperature. Korshunov et al. s argued that short-range 
antiferromagnetic fluctuations are the source of a linear- 
T term in the susceptibility of a two-dimensional Fermi 



liquid which allowed them to obtain good agreement with 
experimental data. 

A very important issue concerning the spectral and 
magnetic properties of the pnictides is the role of 
Coulomb correlations. It is now generally accepted 
that electronic correlations in the pnictides are not as 
strong as in the cuprates and should be classified as 
moderat e) 9 ' 10 . It was show n 10 ' 11 that the spectral proper- 
ties of the pnictides can be reproduced by first-principles 
techniques only if local dynamical Coulomb correlations 
are taken into account. This can be achieved by em- 
ploying the LDA+DMFT method^. This computational 
scheme combines electronic band structure calculations 
in the local density approximation (LDA) with many- 
body physics incorporated in the dynamical mean-field 
theory (DMFT)^. 

In our earlier studj*^ we proposed an explanation of 
the temperature increase of the magnetic susceptibility 
of LaFeAsO which was based on a first-principles analy- 
sis of the low-energy spectral properties caused by local 
dynamical Coulomb correlations, without taking into ac- 
count interatomic magnetic fluctuations. In the present 
paper we employ the LDA+DMFT scheme to demon- 
strate that the proposed mechanism can be applied also 
to understand the origin of the linear temperature de- 
pendence of the uniform magnetic susceptibility in stoi- 
chiometric BaFe2As2- 



II. COMPUTATIONAL METHOD 

The LDA+DMFT computational scheme implemented 
in the present work proceeds in four steps: (i) the 
construction of an effective tight-binding Hamiltonian 
ff WF (k) from a converged LDA solution by projecting 15 
onto Wannier functions, (ii) the addition of the local 
Coulomb interaction H Coul , (iii) a double-counting cor- 
rection which takes into account the local interactions al- 
ready decribed by the LDA, and (iv) the self-consistent 
solution of the DMFT equations on the Matsubara con- 
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tour with continuous-time quantum Monte- Carlo^ (CT- 
QMC) as impurity solver. The Hamiltonian to be solved 
by DMFT is given by 



field, 



H 



DMFT 



(k) = H Wt (k) + H 



rCoul 



H 



DC 



(1) 



Exact investigations of the magnetic response of this 
model must employ a rotationally invariant form of the 
interaction term H Coul . However, up to now there do 
not exist effective algorithms for the solution of a five- 
orbital Hubbard model with the full Coulomb interaction 
within CT-QMC. Furthermore, our investigation of the 
magnetic properties of BaFe2As2 requires a separate, and 
extremely time-consuming, self- consistent DMFT calcu- 
lation for each point of the susceptibility curve. To 
make computations feasible, we therefore include only 
the density-density contributions to the full interaction 
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Here U™> is t ne Coulomb interaction matrix, nf arT is the 
occupation number operator for d electrons in the or- 
bital a or a', with spin a or a', on the ith site. This 
approximation neglects spin flip and pair hopping pro- 
cesses. Nevertheless, as will be shown below, it is able 
to provide correct results for the spectral and magnetic 
properties of BaFe2As2. The double-counting term is 
i? DC = — C/(riDMFT — where 71dmft is the total, self- 
consistently determined number of d electrons obtained 
within LDA+DMFT, and U is the average Coulomb pa- 
rameter for the d shell. 

We construct Wannier functions in the energy window 
including Fe-d and As-p states. Hence, by construction 
energy bands of the H WF Hamiltonian exactly reproduce 
16 Fe-d and As-p bands (two As and Fe atoms in the for- 
mula unit, one formula unit in the unit cell) obtained in 
LDA calculations, and the p—d hybridization is explicitly 
taken into account. 

The interaction matrix U™i is parametrized by the 
effective on-site Coulomb parameter U and intra-orbital 
exchange parameter J according to the procedure de- 
scribed in Ref. 17 . In the present calculation we use 
(7=3.5 eV and 7=0.85 eV obtained with the constrained 
DFT procedure!^. 

The orbitally resolved spectral functions of the inter- 
acting system are then computed as 

A a (u>) = -ilm^[( W +/,)i-^ WF (k)- J ff DC -S( W )]-i, 

k 

(3) 

where the subscript a refers to an orbital, fi is the self- 
consistent chemical potential, and is the self-energy 
on the real axis obtained by analytic continuation using 
the Pade approximan1j2i technique; details are described 
in Refii. 

The uniform magnetic susceptibility is calculated as 
the response of the system to a weak external magnetic 



_ dM{T) d[n t (T)- ni (T)} 



dE h 



dE h 



where M(T) is the field-induced magnetization, n a (T) is 
the number of electrons with spin er, and Eh is the energy 
correction corresponding to the applied field. Since the 
field is finite the calculations are performed in three steps: 
First, we check that the polarization is zero in the absence 
of the field, then we check that M(T) is a linear function 
of E^, and finally we evaluate the derivative in Eq. ^ 
as a ratio of M(T) and E h . 



III. RESULTS 

A. Temperature dependence of the uniform 
magnetic susceptibility 

In Fig. [T] the uniform magnetic susceptibility x(T) 
computed within LDA+DMFT is compared with the ex- 
perimental data of Wang et o7— . In the temperature 
range from 200 K to 600 K the temperature dependence 
of the calculated x(T) is found to be almost perfectly 
linear. However, the slope is by a factor of 1.7 smaller 
than in the experiment. The origin of this quantitative 
discrepancy is not clear at the moment and will be the 
subject of future investigations. To emphasize the linear- 
ity we plot a least-square fit to the last six points of the 
computed data. We note that the observed linear behav- 
ior is, in fact, due to an extended linear region around 
the turning point of %(T) at ss350 K. The obtained x{ T ) 
has a maximum at about 1000 K and decreases for higher 
temperatures. 

More detailed information about the magnetic proper- 
ties of BaFe2As2 can be obtained from an analysis of the 
orbitally resolved contributions Xa(T) (a = xy, yz, xz, 
3z 2 — r 2 , x 2 ~ y 2 ) to the total magnetic susceptibility. In 
Fig. [2] we show the temperature dependence of the Fe 3d 
susceptibilities of BaFe 2 As 2 obtained in LDA+DMFT. 
All contributions have approximately equal slope in the 
temperature interval from 200 K to 500 K. The d xy or- 
bital of Fe provides the largest contribution to the total 
susceptibility. 



B. Connection between magnetic and spectral 
properties 

In our previous paper— we proposed a scenario ac- 
cording to which the anomalous T-behavior of x(T") in 
the pnictides is connected with the presence of a sharp 
peak in the spectral function below the Fermi energy. In 
Fig. [3] the total Fe 3d spectral function A(cj) computed 
within LDA+DMFT is shown in comparison with the 
LDA result. It demonstrates that dynamical correlation 
effects strongly renormalize the spectrum in the vicinity 
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FIG. 1: (Color online) Uniform magnetic susceptibility x{T) 
of BaFe2As2 calculated within LDA+DMFT (squares) in 
comparison with experimental data of Wang et al? (circles 
and triangles). The line is a least-square fit to the last six 
points of the computed data. The inset shows the theoretical 
curve for the full temperature interval. 
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FIG. 2: (Color online) LDA+DMFT results for the orbitally 
resolved Fe 3d susceptibilities Xa(T) with a = xy, yz, xz, 
3z 2 — r 2 , x 2 — y 2 of BaFe2As2 vs. temperature obtained from 
the derivative of the magnetization. 

of the Fermi energy. In particular, the electronic correla- 
tions are seen to lead to a narrow peak below the Fermi 
level while the remaining part of the spectrum is only 
weakly affected by the correlations. The peak in the en- 
ergy window from -4 eV to -2 eV should not be mistaken 
for a lower Hubbard band since a similar peak is already 
present in the LDA result. The Fe 3d spectral weight 
in this energy area is a consequence of the hybridization 
with As 4p states. 

A comparison of the orbitally resolved spectral func- 
tions computed within LDA and LDA+DMFT, respec- 
tively, is shown in Fig.H The LDA+DMFT results are m 
good agreement with previously reported theoretical and 
experimental spectral. Except for the d x 2_ y 2 orbital the 



spectral functions obtained by LDA+DMFT all show a 
sharp peak below the Fermi energy. These peaks origi- 
nate from local correlation effects as pointed out above. 

A quantitative measure of the correlation strength is 
the quasiparticle renormalization factor Z. In the single- 
orbital case it can be obtained from the real-axis self- 
energy S(w) which is related to the effective mass en- 
hancement m*/m by Z~ x = \ — dRe{Yi{uj))/duj = m*/m. 
For a multi-orbital problem the self-energy is a matrix. 
Therefore different orbitals have different m* /m. The 
computed values of the mass enhancement range from 2.5 
to 3.7 and agree well with previous estimates of m* /m for 
pnictides obtained from the renormalization of the LDA 
band structure ^ 1 1 i 22 i 23 . 

To identify possible reasons for the anomalous behavior 
of the magnetic properties it is instructive to compute the 
temperature evolution of the orbitally resolved spectral 
functions. The spectral curves computed for tempera- 
tures ranging from 232 K to 580 K are shown in Fig. [5] 
All spectral functions, except for the d x 2_ y 2 orbital, show 
a temperature sensitive peak which is located approxi- 
mately 100 meV below the Fermi energy. These peaks 
increase in amplitude and become narrow with decreas- 
ing T. We note that the magnitudes of the orbital con- 
tribution \a (T) are proportional to the density of states 
of the corresponding orbitals at the Fermi energy. 

The simplest way to establish a possible connection be- 
tween the temperature evolution of the magnetic suscep- 
tibility and the excitation spectrum in the multi-orbital 
case is to estimate the susceptibility (per spin) using 
the bubble diagram obtained by convoluting the DMFT 
Green's functions, Eq. (JS), 

X°a(T)^^ G aa ,(k,icj)G a , a (k,icj), (5) 
where G aa > (k, iu) = [{iuj + fi)I - F DMFT (k) - ±(iw)] aa > . 
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FIG. 3: (Color online) Comparison of the total Fe 3d spectral 
function of BaFe 2 As 2 as obtained from LDA+DMFT (solid 
line) and LDA (shaded area), respectively. 
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This expression describes the spin susceptibility in the 
absence of vertex corrections, i.e., gives an estimate of 
the magnetic response due to single-particle excitations 
which are characterized by the interacting spectral func- 
tion. An explicit connection between the magnetic re- 
sponse, Eq. ([5]), with the excitation spectrum can be 
made in the one-orbital case when off-diagonal elements 
of G are absent: 



/f(^i) - If {^2) 



A(wi)A(w 2 ). 



( 6 ) 

Here A(w) is the spectral function of the interacting sys- 
tem(per spin), and the temperature enters via the Fermi 
function, /f(oj). In real compounds multi-orbital physics, 
including hybridization effects, is important, in which 
case the full matrix Green functions must be used in 
Eq. ©. 

The influence of interaction effects may be estimated 
by calculating the magnetic susceptibility in the random- 
phase approximation (RPA) . If the orbital dependence of 
the Coulomb interaction between d electrons is neglected, 
the orbital contributions to the total uniform d magnetic 



susceptibility within RPA are given by 

2X°(T) 



X 



RPA 



(T) 



(7) 



where the factor 2 is due to the spin degeneracy. The 
temperature dependence of the orbitally resolved sus- 
ceptibilities computed according to Eq. 0, is shown in 
Fig. [5] For the d 3z 2_ r 2 and d x 2_ y 2 orbitals the results are 
in good qualitative agreement with the full LDA+DMFT 
solution. By contrast, the temperature dependence of the 
susceptibilities corresponding to the d xy , d yzi and d xz 
states is not reproduced by Eq. ([7]); apparently vertex 
corrections are important in this case. 

The above results suggest the following interpretation, 
whose correctness will be demonstrated later: The elec- 
tronic states forming the sharp peak in the spectral func- 
tion below the Fermi energy lead to thermal excitations 
which contribute to the susceptibility. When the energy 
k^T is larger than the distance between the peak and the 
Fermi level, the number of states which can be excited 
is reduced and the susceptibility starts to decrease. A 
more complex mechanism is responsible for the increase 
of the d x 2_ y 2 susceptibility where the corresponding spec- 
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FIG. 4: (Color online) Orbitally resolved Fe 3d spectral func- 
tions of BaFe 2 As 2 obtained within LDA+DMFT (solid lines) 
in comparison with LDA results (filled areas). The Fermi 
energy is set to eV. 
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FIG. 5: (Color online) Fe 3d spectral functions of BaFe2As2 
computed within LDA+DMFT in the temperature range 
232 K - 580 K. The Fermi energy corresponds to eV. 
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FIG. 6: (Color online) LDA+DMFT results for the orbitally 
resolved Fe 3d susceptibilities ~x^ PA (T) with a = xy, yz, xz, 
3z 2 — r 2 , x 2 —y 2 of BaFe2As2 vs. temperature, obtained from 
the RPA expression for the susceptibility in Eq. (Q). 



tral function does not show a peak below the Fermi en- 
ergy. 

Namely, as in the case of LaFeAsOi^ the off- 
diagonal contributions to the susceptibility, Xaa' = 
l//3X)k tu ^aa'(k, iuj)G a i a (k, ioj), and in particular 
y 2 3z 2_ r 2(T), are responsible for the increase of 

y2 (T) . Thus the temperature increase of x^-y 2 (-0 
is caused by the magnetic response of the other orbitals 
(for details see the supplementary material of RefJ^). A 
detailed analysis of the proposed mechanism of the in- 
crease of x(T) using a simplified model will be presented 
in the following section. 
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FIG. 7: (Color online) LDA+DMFT results for the uniform 
magnetic susceptibility of a one-band model with the non- 
interacting Fe d xy DOS of BaFe2 As2 as a function of temper- 
ature. Upper panel: Results computed according to Eq. (j4|; 
the inset shows the spectral function of the interacting system. 
Lower panel: RPA results for the uniform magnetic suscep- 
tibility using Eq. ©; the divergence of the susceptibility for 
/i=0.2 eV at about 300 K is an artifact of the RPA. 



C. Model analysis 

To further clarify the relation between the shape of 
the spectral function and the anomalous temperature be- 
havior of the magnetic properties we will now perform 
a model calculation, where multi-orbital effects are ne- 
glected. As a first step we compute the DMFT spin sus- 
ceptibility for a single-band model, constructed in such 
a way that the non-interacting system has exactly the 
same spectral function as the one for the d xy orbital of 
the tight-binding Hamiltonian i?wF in Eq. ([T]). In the 
upper panel of Fig. [7] the temperature dependence of the 
spin susceptibility computed for the model within DMFT 
is shown for several values of the chemical potential. The 
corresponding spectral functions of the interacting sys- 
tem are presented in the inset of Fig. [7J Depending on 
the peak position there are two characteristically differ- 
ent temperature dependencies of the susceptibility: (i) 
an increase in the low temperature region with a maxi- 
mum at an intermediate temperature followed by a de- 
crease at higher T (/z=0.5 eV, 0.8 eV), and (ii) a mono- 
tonic decrease with temperature (/Lt=0.2 eV). The former 



regime is obtained when the peak is substantially below 
the Fermi energy, the latter regime corresponds to the 
case when the peak is very close to, or right at, the Fermi 
energy. The result of the estimate of the susceptibility 
according to Eq. ([Jj is presented in the lower panel of 
Fig. [Jj It is seen to reproduce all features of the curves 
computed by Eq. (Q}. 

This confirms that the magnetic response of the sys- 
tem is indeed governed by its spectral properties. Thus, 
the driving force of the non-monotonic behavior of the 
susceptibility is the peak below the Fermi energy. In 
particular, the distance between the peak and the Fermi 
level is important. In spite of the fact that the peak in 
the spectral function exists at the level of LDA, its posi- 
tion is too far from the Fermi energy to cause a serious 
increase of the susceptibility. The correlations shift the 
peak towards the Fermi level which results in more pro- 
nounced temperature dependence of x(T). 

To relate the shape of the spectral function with the 
two temperature regimes of the magnetic susceptibility 
even more explicitly, we now compute the temperature 
behavior of X°(^) f° r a system with a density of states 
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FIG. 8: (Color online) Spin susceptibility Xi(T) of a model 
defined by the density of states (|8]) calculated from Eq. © 
for different values of the peak position. The inset shows the 
density of states plotted relative to the Fermi energy (zero 
energy). 

given by a Gauss function with offset from the energy 
axis, 

AH=A + i?exp[(^) 2 ]. (8) 

The parameters in Eq. ((5J were adjusted such that the 
maximum and the width of the function are close to the 
ones obtained with LDA+DMFT for the material spe- 
cific Hamiltonian for BaFe2As2. Results for com- 
puted according to Eq. (|6]) are shown in Fig. [8] The 
temperature behavior of Xi(T) and its evolution upon 
changes of the peak position in the spectral function are 
seen to qualitatively reproduce all features obtained in 
DMFT for the one-band model. This can be viewed as 
a direct indication that the peculiarities observed in the 
anomalous behavior of the spin susceptibility in DMFT 
originate from the shape of the spectral function in the 
vicinity of the Fermi energy. 

Finally we address the question concerning the micro- 
scopic origin of the peaks in the spectral function below 
the Fermi level. In the following analysis we focus on 
the d xy orbital since the peak in its spectral function is 
sharper than those of the other orbitals. The contribu- 
tions of the d xy states to the band structure are shown 
in the left upper panel of Fig. [9] as "fat bands" (i.e., the 
thickness of a band is proportional to the contribution 
of states with selected symmetry). The peak in the d xy 
spectral function (Fig. [SJ right upper panel) is formed 
by regions of relatively flat bands centered at the energy 
approximately -0.4 eV relative to the Fermi level. 

To construct a minimal model describing the energy 
dispersion of the d xy states we solve an effective two- 
band Hamiltonian H 2D (k) with two Fe atoms in the unit 
cell, which is obtained from a projection of Bloch states 
in the vicinity of the Fermi energy onto a subspace of 
Wannier functions with d xy symmetry. The energy bands 
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FIG. 9: (Color online) Band structure and spectral functions 
computed for BaFe2As2; the Fermi energy corresponds to 
eV. (a) Dispersion curves calculated within LDA (dashed 
curves), contributions of the orbitals with d xy symmetry (fat 
bands), energy bands of a two-orbital model obtained as a 
projection onto the d xy states (solid curves), (b) Spectral 
function of the d xy orbital from LDA. (c) Energy bands of 
a model corresponding to the real hopping parameters, (d) 
Spectral function of the model Hamiltonian ([9]). 

of H 2D (k) are shown in the left upper panel of Fig. [9] by 
solid curves. In the next step we introduce the real-space 
Hamiltonian H 2 ® al written for a square lattice with two 
atoms in the unit cell and nonzero hoppings within three 
coordinate spheres. The Hamiltonian H 2 ® al has the form 

Uleai = t E 4^ + f ' E ^R' C * + *" E C lR" ( 9 ) 

where i labels the atoms, and the radius vectors R, R' 
and R" correspond to the cluster of nearest, next-nearest 
and next-next-nearest neighbors, respectively, around 
atom i. The hopping parameters £=-170 meV, t'=98 meV 
and £"=21 meV were computed as Fourier transforms of 
the material-specific Hamiltonian H 2D {k). 

The shape of the energy bands and the spectral func- 
tion computed for Hamiltonian are shown in the lower 
panels of Fig. [§] They are in good agreement with the 
corresponding characteristics obtained from the direct 
calculation. As in the multi-orbital case, the band struc- 
ture of © represents a combination of dispersive bands 
and bands with less pronounced k-dependence. In par- 
ticular, it follows that the peak below the Fermi energy 
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is formed by a relatively flat band located in the same 
energy interval as the peak. 
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FIG. 10: (Color online) Spectral functions of Hamiltonian 
([9]) calculated in the basis of bonding (solid curve) and an- 
tibonding (dashed curve) wave functions. The total spectral 
function is shown by a dotted line. The Fermi energy is set 
to zero. 

In order to understand the coexistence of the flat and 
dispersive regions within a band of given symmetry (d xy 
in our case) it is instructive to plot the spectral functions 
of the Hamiltonian Q in the basis of bonding and an- 
tibonding wave functions. In other words, if is the 
wave function on the one atom in the unit cell, |<fo) is 
the wave function of the other atom. Then the new basis 
is defined as \<f> + ) = \(f>i) + \4> 2 ) and = \4>i) - \(f> 2 ). 
Spectral functions of the Hamiltonian Q computed for 
that basis with realistic hopping parameters and chemical 
potential are shown in Fig. QUI From Fig. [10] it follows 
that the band with a pronounced dispersion is mainly 
formed by antibonding linear combination \<fi—) and the 
main contribution to the peak is provided by the bonding 
function \<f>+) whose dispersion is less pronounced. 



IV. CONCLUSIONS 

By employing the LDA+DMFT method we investi- 
gated the interplay between the spectral and magnetic 
properties of BaFe2 As2 . The calculated temperature de- 
pendence of the uniform magnetic susceptibility is in 
good agreement with experimental data. Our calcula- 
tions show that there are pronounced, temperature sen- 
sitive peaks below the Fermi energy in the spectral func- 
tion of BaFe2As2- We proposed a scenario according to 
which the temperature increase of the susceptibility is 
a consequence of the thermal excitation of the electronic 
states which lead to these peaks. Our analysis is based on 
the DMFT solution of a one-band model with a density 
of states corresponding to the Fe-d xy spectral function of 
BaFe2As2- The results clearly demonstrate that the peak 
in the spectral function in the vicinity of the Fermi en- 
ergy is a prerequisite for the linear temperature increase 
of magnetic susceptibility. The peaks in the real com- 
pound are due to the weak dispersion of bonding states 
arising from the layered structure of BaFe2 As2 . 
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